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Abstract 

We evaluate the dipolar oscillations of a harmonically trapped fermion gas contain- 
ing thermal bosonic impurities as a function of the anisotropy of the trap, from the 
numerical solution of the Vlasov-Landau equations for the one-body phase-space 
distribution functions. Starting from a situation in which the two components of 
the gaseous mixture perform almost independent oscillations inside a spherical trap, 
we demonstrate that different collision behaviors arise for oscillations in different 
directions as the trap is deformed into an elongated cigar-like shape. An increase 
in the anisotropy of the confinement thus suffices to drive a transition of dipolar 
modes from a collisionless to a collisional regime. 
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1 Introduction 



Many aspects of the physics of ultracold quantum gases and their quantum 
phase transitions can usefully be explored by measuring the frequencies and 
the damping rates of monopolar "breathing" modes and of quadrupolar modes 
at various temperatures (see Ref. [1] and references therein). These measure- 
ments, which have been performed with rather high accuracy on both bosonic 
and fermionic systems, have been important to the understanding of Bose- 
Einstein condensation and of degenerate and superfluid fermion gases in con- 
fined geometries, and have become a testing ground for many-body theories. 

In contrast, the dipolar "sloshing" modes in a monatomic gas under harmonic 
confinement contain no deep information on its microscopic behavior. These 
modes obey a theorem stemming from the work of Kohn [5] on a fiuid of 
charged particles with arbitrary momentum-conserving interactions in a mag- 
netic field, which was generalized by later workers [6] by adding an external 
scalar harmonic potential extending throughout all space. The theorem states 
that in a system of interacting particles confined by an external harmonic po- 
tential T4(r) = r ■ K • r/2 the dynamics of the center of mass is completely 
decoupled from that of the internal degrees of freedom. This implies sharp 
resonances at the bare trap frequencies for a uniform exciting field, which 
can only rigidly drive the density distribution of the system. The so-called 
extended Kohn theorem (EKT) has had important consequences in the de- 
velopment of time-dependent density functional theory for the dynamics of 
electrons in confined geometries such as quantum wells and quantum dots [7] . 

Turning, however, to a gaseous mixture, it is evident that its dipolar modes 
provide a direct signature of its coUisionality [8,9]. In a regime of low collision 
rates the components of the mixture are independently oscillating, whereas 
they oscillate together when their coUisionality is sufficiently high. This dy- 
namical transition from a collisionless to a collisional (hydrodynamic) behav- 
ior has been followed both experimentally and numerically in two-component 
fermion mixtures [8,9]. The role played in this context by mobile impurities 
inside a fermion gas under spherical confinement has also been studied by 
numerical means [10]. Interatomic collisions also play an important role in 
the development and understanding of techniques for the cooling of gases in 
essentially harmonic confinements down to ultralow temperatures [2,3,4]. 

In the present work we examine the role of the anisotropy of the harmonic 
confinement in this dynamical transition for a fermion-boson mixture simu- 
lating a gas of fermionic '^'^K atoms which contains a small concentration of 
thermal ^''Rb atoms. We numerically solve the Vlasov-Landau equations for 
the evolution of the phase-space distribution functions within a particle-in-cell 
approach (see Refs. [10,11] for the technical details). In particular, our results 
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demonstrate that in a cigar-shaped harmonic confinement, rigid oscillations 
of the centre of mass of the whole mixture in the "soft" axial direction can 
coexist with independent oscillations of the two components in a "hard" radial 
direction. This result can be considered as a further extension of the EKT. 

The paper is organized as follows. In Sec. 2 we describe the mixture under 
study, while in Sees. 3 and 4 we analyze its dynamics under time-dependent 
drives. Finally, Sec. 5 offers some concluding remarks. 



2 The physical system 

The system is a mixture of spin-polarized atoms and ^^Rb atoms inside 
a magnetic trap and the interspecies interaction is modeled by a strongly 
repulsive contact potential characterized by an s-wave scattering length a = 
2 X 10'^ Bohr radii. The mixture is at low temperature and the number of ^''Rb 
bosons is kept low in order to avoid the formation of a condensate, which may 
in turn reduce the coUisionality of the system [12]. Moreover, this allows us to 
neglect the boson-boson interactions. 

Due to the difference in masses and in internal states, fermions and bosons 
experience different confining potentials. These are given by 

^F,B(r) = ^ mF,B c^F,B(r' + A' z^) , (1) 

where we have assumed the same anisotropy parameter A for the two species. 
In Eq. (1) mi;' and rriB are the atomic masses of fermions and bosons and uup 
and cub are their trap frequency along the radial direction. The effects of the 
anisotropy will be examined by varying A from 1 to 0.01 while keeping fixed the 
average cDf.b = ujf,b^^^^ of the trap frequencies. This value will be taken as the 
geometric average of the frequencies in the experiments carried out at LENS 
on the same mixture [13]. In this way we intend to explore the dependence 
of the dynamics solely on the anisotropy of the traps without changing other 
properties of the gas such as the mean densities and the chemical potentials. 
For definiteness we also set the temperature of the sample at T = 0.2Tf with 
Tp = fvjjFiQNpY^^ /ks being the Fermi temperature and Np the number of 
fermions. 



3 ColUsionality near equiUbrium 

The number of collisions per unit time {i.e. the collision rate) in the mixture 
at equilibrium is both classically and quantum-mechanically independent of A 
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for fixed Oi>F,B- However, a slight mismatch of the trap centers can lead to a 
strong decrease in the collisionality, which can be understood to be the result 
of the decrease in the overlap of the two density profiles. In the classical limit, 
the coUisionahty in a mixture of Np fermions and Nb bosons reads 



when the trap centers have been offset by {xq, zq). Here P = I/ZcbT, a — Atto^ 
is the total scattering cross section, rrir the reduced boson-fermion mass, Kr = 
KpKb/ {Kp + Kb) the effective oscillator strength, and Kp^B = ^p^bi^f^b 
the average oscillator strength of each species. Taking xq and Zq to be of 
comparable magnitudes, the collisionality as a function of the anisotropy thus 
displays two different trends depending on the direction of the offset. 

In the numerical simulations we first shift rigidly the fermionic cloud in ei- 
ther the radial or the axial direction by Uho — {h/mpupy^^, which is usually 
much smaller than the radius of the two clouds. Then the fermion cloud starts 
oscillating around its equilibrium position and collisions among bosons and 
fermions can occur. During the evolution the relative position of the centers 
of the two clouds changes continuously and Eq. (2) gives only an indication 
of the dependence of the time-averaged coUisionality on A. Prom the simula- 
tion data we have calculated the average number of collisions per unit time 
during a fixed time interval, that we choose to be t = lOn/ujp. The result is 
shown in Fig. 1, where we plot the average quantum collision rates and 
Fg, corresponding to dipolar oscillations of the fermions in the two directions, 
as functions of the trap anisotropy A. The filled symbols refer to the collision 
rates in units of the constant frequency scale Cup. 

The data in Fig. 1 show that for the system parameters that we are consid- 
ering the variation with A of the two collision rates F^ in these absolute units 
becomes appreciable only at low A. However, referring to the empty symbols in 
Fig. 1 which report F^ and F^ in units of the fermionic axial and radial frequen- 
cies = X^^^ujp and Ur = \~^^^ujp, it is seen that the numbers of collisions 
made on average during an oscillation period in the axial or in the radial di- 
rection are becoming very different. This suggests that, since Tg/cOr ^ 1 while 
I'q/ujz is strongly increasing, the mixture may have a coUisionless behavior 
along the radius r and a hydrodynamic behavior along the z axis. 

The return of the system to a new equilibrium state (including instantaneous 
mean-field effects) can be directly investigated in the Vlasov-Landau frame by 
following the collision rate as a function of time. As the evolution proceeds 
the average number of collisions approaches that corresponding to the new 
equilibrium. We display in Fig. 2 the time dependence of the radial collision 
rate for different values of A. The collision rate diminishes relatively rapidly 
towards the final equilibrium state, as can be expected from the temperature 
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Fig. 1. Total collision rates Fg per fermion in the radial (i = r) and axial (i = z) 
directions, as functions of the anisotropy parameter A. Filled symbols show the 
collision rates in units of top and empty symbols report V^^/ {NpuJi). The squares 
refer to oscillations along the radial direction and the circles to oscillations along 
the axis of the trap. The lines are guides to the eyes. 
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Fig. 2. Total collision rate Fg(t)/rg(0) as a function of time for several values of 
the anisotropy A and in the case of an initial radial displacement. 



dependence of T^x in Eq. (2) since the system performs a number-conserving 
evolution. 
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4 Collisionless vs collisional behaviour 



Let us summarize some general features of the dynamical behavior of the mix- 
ture that are associated with an initial displacement of the fermion cloud, 
before turning to discuss the role of trap anisotropy. In the absence of inter- 
actions the fermions oscillate without damping at their bare trap frequency 
while the bosonic impurities remain at rest. On switching on the interactions 
the fermionic components sets the bosons into motion and beats between os- 
cillations at the bosonic and the fermionic trap frequencies appear [11]. If the 
collision rate is much lower than the trap frequencies, the damping rate is 
small and the two clouds oscillate for many periods. On increasing the col- 
lisionality up to a regime where the collision rate is of the same order as 
the trap frequencies, the oscillation frequencies are not well defined and the 
center-of-mass motions of the two clouds become strongly damped. Finally, 
at very large coUisionality the two clouds become glued together and oscillate 
at the same frequency without noticeable damping. Thus a clear signature of 
the onset of a collisional regime is the locking of the oscillations of the two 
components and the vanishing of their damping rates. 

We turn now to examining the role played by trap anisotropy in this scenario. 
We plot in Fig. 3(a) the oscillation frequencies of the fermions in the radial 
and axial direction as functions of the anisotropy parameter A which drives 
the collision rates referred to the trap frequencies (see Fig. 1). The data in 
Fig. 3 have been obtained by fitting the oscillatory motions of the fermions 
with functions of the form cos(f)t-|-0) exp(— 7t). The oscillation frequency fir 
in the radial direction approaches the radial trap frequency ujr for low A and 
decreases below this value with increasing A. On the contrary, the oscillation 
frequency Qz iii the axial direction takes its maximum value at A = 1 and 
decreases with decreasing A. These two opposite trends can be understood if 
we plot the oscillation frequencies as functions of the collision rates scaled by 
the trap frequencies. In such a plot (see Fig. 3(b)) both oscillation frequencies 
lie on a single monotonically decreasing curve. The plot suggests that the 
leftmost data correspond to collisionless radial oscillations while the rightmost 
data reflect an approach to the collisional regime where the axial oscillation 
frequencies are close to the hydrodynamic value. 

To support the idea that the mixture at large trap anisotropics is in different 
regimes in the radial and in the axial direction, we report in Fig. 4 the damping 
rates 7^ and ^y^ for the oscillatory motions as functions of the corresponding 
collision rates. In the radial direction (circles in Fig. 4) the damping rate 
increases as a function of the scaled collision rate F^/cj^ and seems to reach a 
maximum: this is a signature that the gas is moving from the collisionless to the 
intermediate collisional regime with decreasing anisotropy. On the contrary, 
in the axial direction (squares in Fig. 4) the damping rate decreases as a 
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Fig. 3. Radial and axial oscillation frequencies (in units of oor and a;^, respectively) 
as functions of (a) the anisotropy parameter A, and (b) the collision rate V^/Np per 
fermion (in units of uJr for radial oscillations and of Uz for axial oscillations). 
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Fig. 4. Radial (circles) and axial (squares) damping rates 7j (in units of cOr and w^, 
respectively) as functions of the collision rates T^g/{Npu;r) and T^/{NpuJz)- 

function ofT^/uz and this is a typical behavior of an intermediate-to-collisional 
transition. 



5 Conclusions 



We have studied the dipolar modes of a harmonically confined fermion-boson 
mixture as functions of the trap anisotropy and found that, although the 
collisionality is independent of the anisotropy, a collisionless behavior in the 
radial direction and a hydrodynamic behavior in the axial direction can be 
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simultaneously established for large anisotropies. This conclusion is based on 
the analysis of the dipolar oscillations of a fermion gas at ultralow temperature 
and containing a limited concentration of thermal bosons. 

The analysis that we have presented has only regarded the dipolar oscillations 
of the gas. Of course, other oscillatory modes and the dynamics of a ballistic 
expansion would provide further insight on the physical behavior of the gas. 
The results of these further studies will be presented elsewhere. 
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